Genome-wide association study identifies DRAM1 associated with Tourette syndrome in Taiwan.

BACKGROUND: Tourette syndrome (TS) is a neurodevelopmental disorder characterized by motor and vocal tics. Several susceptibility loci associated with TS have been identified previously in populations of European descent using genome-wide association studies (GWAS). However, the exact pathogenic mechanism underlying TS is unknown; additionally, the results of previous GWAS for TS were based on Western populations, which may not translate to other populations. Therefore, we conducted a GWAS in Taiwanese patients with TS and chronic tic disorders (CTDs), with an aim to elucidate the genetic basis and potential risk factors for TS in this population. METHODS: GWAS was performed on a Taiwanese TS/CTDs cohort with a sample size of 1,007 patients with TS and 25,522 ancestry-matched controls. Additionally, polygenic risk score was calculated and assessed. RESULTS: Genome-wide significant locus, rs12313062 (p=1.43 × 10-8) and other 9 single nucleotide polymorphisms, were identified in chromosomes 12q23.2, associated with DRAM1 and was a novel susceptibility locus identified in TS/CTDs group. DRAM1, a lysosomal transmembrane protein regulated by p53, modulates autophagy and apoptosis, with potential implications for neuropsychiatric conditions associated with autophagy disruption. CONCLUSIONS: This study conducted the first GWAS for TS in a Taiwanese population, identifying a significant locus on chromosome 12q23.2 associated with DRAM1. These findings provide novel insights into the neurobiology of TS and potential directions for future research in this area.

Utility of polygenic scores across diverse diseases in a hospital cohort for predictive modeling.

Polygenic scores estimate genetic susceptibility to diseases. We systematically calculated polygenic scores across 457 phenotypes using genotyping array data from China Medical University Hospital. Logistic regression models assessed polygenic scores' ability to predict disease traits. The polygenic score model with the highest accuracy, based on maximal area under the receiver operating characteristic curve (AUC), is provided on the GeneAnaBase website of the hospital. Our findings indicate 49 phenotypes with AUC greater than 0.6, predominantly linked to endocrine and metabolic diseases. Notably, hyperplasia of the prostate exhibited the highest disease prediction ability (P value = 1.01 × 10-19, AUC = 0.874), highlighting the potential of these polygenic scores in preventive medicine and diagnosis. This study offers a comprehensive evaluation of polygenic scores performance across diverse human traits, identifying promising applications for precision medicine and personalized healthcare, thereby inspiring further research and development in this field.

Menarche-a journey into womanhood: age at menarche and health-related outcomes in East Asians.

STUDY QUESTION: Are there associations of age at menarche (AAM) with health-related outcomes in East Asians? SUMMARY ANSWER: AAM is associated with osteoporosis, Type 2 diabetes (T2D), glaucoma, and uterine fibroids, as demonstrated through observational studies, polygenic risk scores, genetic correlations, and Mendelian randomization (MR), with additional findings indicating a causal effect of BMI and T2D on earlier AAM. WHAT IS KNOWN ALREADY: Puberty timing is linked to adult disease risk, but research predominantly focuses on European populations, with limited studies in other groups. STUDY DESIGN, SIZE, DURATION: We performed an AAM genome-wide association study (GWAS) with 57 890 Han Taiwanese females and examined the association between AAM and 154 disease outcomes using the Taiwanese database. Additionally, we examined genetic correlations between AAM and 113 diseases and 67 phenotypes using Japanese GWAS summary statistics. PARTICIPANTS/MATERIALS, SETTING, METHODS: We performed AAM GWAS and gene-based GWAS studies to obtain summary statistics and identify potential AAM-related genes. We applied phenotype, polygenic risk scores, and genetic correlation analyses of AAM to explore health-related outcomes, using multivariate regression and linkage disequilibrium score regression analyses. We also explored potential bidirectional causal relationships between AAM and related outcomes through univariable and multivariable MR analyses. MAIN RESULTS AND THE ROLE OF CHANCE: Fifteen lead single-nucleotide polymorphisms and 24 distinct genes were associated with AAM in Taiwan. AAM was genetically associated with later menarche and menopause, greater height, increased osteoporosis risk, but lower BMI, and reduced risks of T2D, glaucoma, and uterine fibroids in East Asians. Bidirectional MR analyses indicated that higher BMI/T2D causally leads to earlier AAM. LIMITATIONS, REASONS FOR CAUTION: Our findings were specific to Han Taiwanese individuals, with genetic correlation analyses conducted in East Asians. Further research in other ethnic groups is necessary. WIDER IMPLICATIONS OF THE FINDINGS: Our study provides insights into the genetic architecture of AAM and its health-related outcomes in East Asians, highlighting causal links between BMI/T2D and earlier AAM, which may suggest potential prevention strategies for early puberty. STUDY FUNDING/COMPETING INTEREST(S): The work was supported by China Medical University, Taiwan (CMU110-S-17, CMU110-S-24, CMU110-MF-49, CMU111-SR-158, CMU111-MF-105, CMU111-MF-21, CMU111-S-35, CMU112-SR-30, and CMU112-MF-101), the China Medical University Hospital, Taiwan (DMR-111-062, DMR-111-153, DMR-112-042, DMR-113-038, and DMR-113-103), and the Ministry of Science and Technology, Taiwan (MOST 111-2314-B-039-063-MY3, MOST 111-2314-B-039-064-MY3, MOST 111-2410-H-039-002-MY3, and NSTC 112-2813-C-039-036-B). The funders had no influence on the data collection, analyses, or conclusions of the study. No conflict of interests to declare. TRIAL REGISTRATION NUMBER: N/A.

Identification of Genetic Variants for Risk Prediction and Early Diagnosis of Age-Related Macular Degeneration in the Taiwanese Population.

Age-related macular degeneration (AMD) is the leading cause of blindness in the elderly worldwide. The prevalence and phenotypes of AMD differ among populations, including between people in Taiwan and other regions. We performed a genome-wide association study to identify genetic variants and to develop genetic models to predict the risk of AMD development and progression in the Taiwanese population. In total, 4039 patients with AMD and 16,488 non-AMD controls (aged ≥ 65 years) were included. We identified 31 AMD-associated variants (p < 5 × 10-8) on chromosome 10q26, surrounding PLEKHA1-ARMS2-HTRA1. Two genetic models were constructed using the clump and threshold method. Model 1 included the single nucleotide polymorphism rs11200630 and showed a 1.31-fold increase in the risk of AMD per risk allele (95% confidence interval (CI) = 1.20-1.43, p < 0.001). In model 2, 1412 single-nucleotide polymorphisms were selected to construct a polygenic risk score (PRS). Individuals with the top 5% PRS had a 1.40-fold higher AMD risk compared with that of individuals with a PRS in the bottom quartile (95% CI = 1.04-1.89, p = 0.025). Moreover, the PRS in the upper quartile was related to a decreased age at AMD diagnosis by 0.62 years (95% CI = -1.15, -0.09, p = 0.023). Both genetic models provide useful predictive power for populations at high risk of AMD, affording a basis for identifying patients requiring close follow-up and early intervention.

Genetic susceptibility to prostate cancer in Taiwan: A genome-wide association study.

We conducted the first genome-wide association study (GWAS) of prostate cancer (PCa) in Taiwan with 1844 cases and 80,709 controls. Thirteen independent single-nucleotide polymorphisms (SNPs) reached genome-wide significance (p < 5 × 10-8 ). Among these, three were distinct from previously identified loci: rs76072851 in CORO2B gene (15q23), odds ratio (OR) = 1.54, 95% confidence interval (CI), 1.36-1.76, p = 5.30 × 10-11 ; rs7837051, near two long noncoding RNA (lncRNA) genes, PRNCR1 and PCAT2 (8q24.21), OR = 1.41 (95% CI, 1.31-1.51), p = 8.77 × 10-21 ; and rs56339048, near an lncRNA gene, CASC8 (8q24.21), OR = 1.25 (95% CI, 1.16-1.35), p = 2.14 × 10-8 . We refined the lead SNPs for two previously identified SNPs in Taiwanese: rs13255059 (near CASC8), p = 9.02 × 10-43 , and rs1456315 (inside PRNCR1), p = 4.33 × 10-42 . We confirmed 35 out of 49 GWAS-identified East Asian PCa susceptibility SNPs. In addition, we identified two SNPs more specific to Taiwanese than East Asians: rs34295433 in LAMC1 (1q25.3) and rs6853490 in PDLIM5 (4q22.3). A weighted genetic risk score (GRS) was developed using the 40 validated SNPs and the area under the receiver-operating characteristic curve for the GRS to predict PCa was 0.67 (95% CI, 0.63-0.71). These identified SNPs provide valuable insights into the molecular mechanisms of prostate carcinogenesis in Taiwan and underscore the significant role of genetic susceptibility in regional differences in PCa incidence.

Identification of the efficacy of parentage testing based on bi-allelic autosomal single nucleotide polymorphism markers in Taiwanese population.

Parentage testing is crucial for forensic DNA analysis, using short tandem repeats (STRs). Single nucleotide polymorphisms (SNPs) with high minor allele frequency (MAF) are promising for human identification. This study aimed to develop SNP markers for parentage testing in the Taiwanese population and compare their accuracy with STRs. The TPMv1 SNP microarray (714,457 SNPs) was used to screen 180,000 Taiwanese individuals and analyze the SNP data using PLINK. After quality control, allelic distribution, and MAF considerations, a set of SNPs with significant inheritance information was selected. Parentage testing was conducted on 355 single parent-child pairs using both STRs and SNPs, employing three kinship algorithms: identity by descent, kinship-based inference for genome-wide association studies, and the combined paternity index/probability of paternity (CPI/PP). An Affymetrix signature probe for kinship testing (ASP) was also used. Based on the quality control and selection criteria, 176 SNPs with MAF > 0.4995 were selected from the Taiwanese population. The CPI/PP results calculated using SNPs were consistent with the STR results. The accuracy of the SNPs used in the single-parent-child parentage testing was > 99.99%. The set of 176 SNPs had a higher identification rate in the single parent-child parentage test than in the ASP. The CPI/PP value calculated using 176 SNPs was also more accurate than that calculated using ASP. Our findings suggest that these 176 SNPs could be used for single-parent-child parentage identification in the Taiwanese population.

The genome-wide association study of serum IgE levels demonstrated a shared genetic background in allergic diseases.

Immunoglobulin E (IgE) synthessis is highly related to a variety of atopic diseases, and several genome-wide association studies (GWASs) have demonstrated the association between genes and IgE level. In this study, we conducted the largest genome-wide association study of IgE involving a Taiwanese Han population. Eight independent variants exhibited genome-wide significance. Among them, an intronic SNP of CD28, rs1181388, and an intergenic SNP, rs1002957030, on 11q23.2 were identified as novel signals for IgE. Seven of the loci were replicated successfully in a meta-analysis using data on Japanese population. Among all the human leukocyte antigen (HLA) regions, HLA-DQA1*03:02 - HLA-DQB1*03:03 was the most significant haplotype (OR = 1.25, SE = 0.02, FDR = 1.6 × 10-14), corresponding to HLA-DQA1 Asp160 and HLA-DQB1 Leu87 amino acid residues. The genetic correlation showed significance between IgE and allergic diseases including asthma, atopic dermatitis, and pollinosis. IgE PRS was significantly correlated with total IgE levels. Furthermore, the top decile IgE polygenic risk score (PRS) group had the highest risk of asthma for the Taiwan Biobank and Biobank Japan cohorts. IgE PRS may be used to aid in predicting the occurrence of allergic reactions before symptoms occur and biomarkers are detectable. Our study provided a more comprehensive understanding of the impact of genomic variants, including complex HLA alleles, on serum IgE levels.

The largest genome-wide association study for breast cancer in Taiwanese Han population.

PURPOSE: Breast cancer is a molecularly heterogeneous disease, and multiple genetic variants contribute to its development and prognosis. Most of previous genome-wide association studies (GWASs) and polygenic risk scores (PRSs) analyses focused on studying breast cancers of Caucasian populations, which may not be applicable to other population. Therefore, we conducted the largest breast cancer cohort of Taiwanese population to fill in the knowledge gap. METHODS: A total of 152,534 Participants recruited by China Medical University Hospital between 2003 and 2019 were filtered by several patient selection criteria and GWAS quality control steps, resulting in the inclusion of 2496 cases and 9984 controls for this study. We then conducted GWAS for all breast cancers and PRS analyses for all breast cancers and the four breast cancer subtypes, including luminal A, luminal B, basal-like, and HER2-enriched. RESULTS: The GWAS analyses identified 113 SNPs, 50 of which were novel. The PRS models for all breast cancers and the luminal A subtype showed positively correlated trends between the PRS and the risk of developing breast cancer. The odds ratios (95% confidence intervals) for the groups with the highest PRS in all breast cancers and the luminal A subtype were 5.33 (3.79-7.66) and 3.55 (2.13-6.14), respectively. CONCLUSION: In summary, we explored the association of genetic variants with breast cancer in the largest Taiwanese cohort and developed two PRS models that can predict the risk of developing any breast cancer and the luminal A subtype in Taiwanese women.

Identification of 13 Novel Loci in a Genome-Wide Association Study on Taiwanese with Hepatocellular Carcinoma.

Liver cancer is caused by complex interactions among genetic factors, viral infection, alcohol abuse, and metabolic diseases. We conducted a genome-wide association study and polygenic risk score (PRS) model in Taiwan, employing a nonspecific etiology approach, to identify genetic risk factors for hepatocellular carcinoma (HCC). Our analysis of 2836 HCC cases and 134,549 controls revealed 13 novel associated loci such as the FAM66C gene, noncoding genes, liver-fibrosis-related genes, metabolism-related genes, and HCC-related pathway genes. We incorporated the results from the UK Biobank and Japanese database into our study for meta-analysis to validate our findings. We also identified specific subtypes of the major histocompatibility complex that influence both viral infection and HCC progression. Using this data, we developed a PRS to predict HCC risk in the general population, patients with HCC, and HCC-affected families. The PRS demonstrated higher risk scores in families with multiple HCCs and other cancer cases. This study presents a novel approach to HCC risk analysis, identifies seven new genes associated with HCC development, and introduces a reproducible PRS model for risk assessment.

Genome-Wide Association Study of Alopecia Areata in Taiwan: The Conflict Between Individuals and Hair Follicles.

Purpose: Alopecia areata (AA) is one of the most prevalent autoimmune diseases affecting humans. Given that hair follicles are immune-privileged, autoimmunity can result in disfiguring hair loss. However, the genetic basis for AA in the Taiwanese population remains unknown. Materials and Methods: A genome-wide association study was conducted using a cohort of 408 AA cases and 8167 controls. To link variants to gene relationships, we used 882 SNPs (P<1E-05) within 74 genes that were associated with AA group to build the biological networks by IPA software. HLA diplotypes and haplotypes were analyzed using Attribute Bagging (HIBAG)-R package and chi-square analysis. Results: Seven single nucleotide polymorphisms (SNPs) including LINC02006 (rs531166736, rs187306735), APC (rs112800832_C_CAT), SRP19 (rs139948960, rs144784670), EGFLAM (rs16903975) and LDLRAD3 (rs79874564) were closely associated with the AA phenotype (P<5E-08). Examination of biological networks revealed that these genomic areas are associated with antigen presentation signaling, B cell and T cell development, Th1 and Th2 activation pathways, Notch signaling, crosstalk signaling between dendritic cells and natural killer cells, and phagosome maturation. Based on human leukocyte antigen (HLA) genotype analysis, four HLA genotypes (HLA-B*15:01-*40:01, HLA-DQA1*01:02-*03:03, HLA-DQA1*01:02, and HLA-DQB1*02:01) were found to be associated with AA (adjusted p-value<0.05). HLA-DQA1*01:02 is the most significantly related gene in the Taiwanese population (adjusted p-value = 2.09E-05). Conclusion: This study successfully identified susceptibility loci associated with AA in the Taiwanese population. These findings not only shed light on the origins of AA within the Taiwanese context but also contribute to a comprehensive understanding of the genetic factors influencing AA susceptibility.

Effects of Human Leukocyte Antigen DRB1 Genetic Polymorphism on Anti-Cyclic Citrullinated Peptide (ANTI-CCP) and Rheumatoid Factor (RF) Expression in Rheumatoid Arthritis (RA) Patients.

Rheumatoid arthritis (RA) is a systemic disease characterized by non-infectious inflammation of the joints and surrounding tissues, which can cause severe health problems, affect the patient's daily life, and even cause death. RA can be clinically diagnosed by the occurrence of blood serological markers, rheumatoid factor (RF) and anti-cyclic citrullinated peptide antibody (anti-CCP). However, about 20% of RA patients exhibit negative results for both markers, which makes RA diagnosis difficult and, therefore, may delay the effective treatment. Previous studies found some evidence that human leukocyte antigen (HLA)-related genes might be the susceptibility genes for RA and their polymorphisms might contribute to varieties of susceptibility and disease severity. This study aimed for the genetic polymorphisms of the RA patient genome and their effects on the RA patient's serological makers, RF and anti-CCP. A total of 4580 patients' electronic medical records from 1992 to 2020 were retrieved from the China Medical University Hospital database. The most representative single-nucleotide polymorphisms (SNPs) were identified through a genome-wide association study (GWAS) followed by enzyme-linked immunosorbent assay (ELISA) validation using the blood from 30 additional RA patients. The results showed significant changes at the position of chromosome 6 with rs9270481 being the most significant locus, which indicated the location of the HLA-DRB1 gene. Further, patients with the CC genotype at this locus were more likely to exhibit negative results for RF and anti-CCP than those with the TT genotype. The C allele was also more likely to be associated with negative results for RF and anti-CCP. The results demonstrated that a genetic polymorphism at rs9270481 affected the expression of RF and anti-CCP in RA patients, which might indicate the necessity to develop a personalized treatment plan for each individual patient based on the genetic profile.

Combining polygenic risk scores and human leukocyte antigen variants for personalized risk assessment of type 1 diabetes in the Taiwanese population.

AIMS: To analyse the genome-wide association study (GWAS) data of patients with type 1 diabetes mellitus (T1D) in order to develop a risk score for the genetic effects on T1D risk and age at diagnosis in the Taiwanese population. MATERIALS AND METHODS: We selected 610 patients with T1D and 2511 healthy individuals from an electronic medical record database of more than 300 000 individuals with genetic information, analysed their GWAS data, and developed a polygenic risk score (PRS). RESULTS: The PRS, based on 149 selected single-nucleotide polymorphisms, could effectively predict T1D risk. A PRS increase was associated with increased T1D risk (odds ratio [OR] 2.09, 95% confidence interval [CI] 1.72-2.55). Moreover, a 1-unit increase in standardized T1D PRS decreased the age at diagnosis by 0.74 years. Combined PRS and human leukocyte antigen (HLA) DQA1*03:02-DQA1*05:01 genotypes could accurately predict T1D risk. In multivariable models, HLA variants and PRS were independent risk factors for T1D risk (OR 3.76 [95% CI 1.54-9.16] and 1.71 [95% CI 1.37-2.13] for HLA DQA1*03:02-DQA1*05:01 and PRS, respectively). In a limited study population of those aged ≤18 years, PRS remained significantly associated with T1D risk. The association between T1D PRS and age at diagnosis was more obvious among males and patients aged ≤18 years. CONCLUSIONS: Polygenic risk score and HLA variations enable personalized risk estimates, enhance newborn screening efficiency for ketoacidosis prevention, and addresses the gap in data on T1D prediction in isolated Asian populations.

Efficacy of Warfarin Therapy Guided by Pharmacogenetics: A Real-world Investigation Among Han Taiwanese.

PURPOSE: The anticoagulation activity of warfarin in populations with CYP2C9, VKORC1, and CYP4F2 variants differs between individuals and is correlated with poor international normalized ratio (INR) control. Pharmacogenetics-guided warfarin dosing has been successfully developed for patients with genetic variations in recent years. However, few real-world data have been used to investigate the INR and warfarin dosage and the time to target INR. This study examined the largest collection of genetic and clinical real-world data related to warfarin to provide further evidence supporting the benefits of pharmacogenetics in clinical outcomes. METHODS: We retrieved a total of 69,610 INR-warfarin records after the index date from 2,613 patients in the China Medical University Hospital database between January 2003 and December 2019. Each INR reading was obtained from the latest laboratory data after the hospital visit date. Patients with a history of malignant neoplasms or pregnancy before the index date were excluded, as were patients without data on INR measurements after the fifth day of prescription, genetic information, or gender variables. The primary outcomes were the INR and warfarin dosage during days 7, 14, 28, 56, and 84 after prescription. The secondary outcome was the time required to reach the INR ranges of 1.5 to 3.0 and >4.0. FINDINGS: A total of 59,643 INR-warfarin records from 2188 patients were retrieved. The average INR was higher for homozygous carriers of the minor allele at CYP2C9 and VKORC1 during the first 7 days (1.83 [1.03] [CYP2C9*1] and 2.46 [1.44] [CYP2C9*3], P < 0.001; 1.39 [0.36] [rs9923231 G/G], 1.55 [0.79] [rs9923231 G/A], and 1.96 [1.13] [rs9923231 A/A], P < 0.001) than for the wild-type allele. These patients with variants required lower warfarin doses than those with the wild-type allele during the first 28 days. CYP4F2 variant patients seemed to require higher doses of warfarin than those in the wild-type group; however, no significant difference in the average INR was observed (1.95 [1.14] [homozygous V433 carriers], 1.78 [0.98] [heterozygous V433M carriers], and 1.66 [0.91] [homozygous M433 carriers], P = 0.016). IMPLICATIONS: Our study indicates that genetic variants in the Han population may enhance warfarin responsiveness, which holds clinical relevance. An increased warfarin dosage was not linked to a shorter time to therapeutic INR between CYP4F2 variant patients and those with a wild-type allele. Assessing CYP2C9 and VKORC1 genetic polymorphisms before initiating warfarin treatment in real-world practice is essential for potentially vulnerable patients and is likely to optimize therapeutic dosing.

Liraglutide With Metformin Therapy Ameliorates Hepatic Steatosis and Liver Injury in a Mouse Model of Non-alcoholic Steatohepatitis.

BACKGROUND/AIM: Non-alcoholic fatty liver disease is a major cause of liver-related morbidity and mortality. Metformin is a widely used medication and may have additional benefits beyond glycemic control. Liraglutide, a novel treatment for diabetes and obesity, also has beneficial effects on non-alcoholic steatohepatitis (NASH). Metformin and liraglutide have both benefited NASH treatment. However, no study has reported the effects of combination therapy with liraglutide and metformin on NASH. MATERIALS AND METHODS: We investigated the in vivo effects of metformin and liraglutide on NASH in a methionine/choline-deficient (MCD) diet-fed C57BL/6JNarl mouse model. Serum triglyceride, alanine aminotransferase and alanine aminotransferase levels were documented. Histological analysis was performed according to the NASH activity grade. RESULTS: After treatment with liraglutide and metformin, body weight loss improved, and the liver/body weight ratio decreased. The metabolic effects and liver injury improved. Liraglutide and metformin alleviated MCD-induced hepatic steatosis and injury. Histological analysis revealed that NASH activity was reduced. CONCLUSION: Our results provide evidence for the anti-NASH activity of liraglutide in combination with metformin. Liraglutide with metformin may offer the potential for a disease-modifying intervention for NASH.

A Genome-Wide Association Study Identified Novel Genetic Susceptibility Loci for Oral Cancer in Taiwan.

Taiwan has the highest incidence rate of oral cancer in the world. Although oral cancer is mostly an environmentally induced cancer, genetic factors also play an important role in its etiology. Genome-wide association studies (GWAS) have identified nine susceptibility regions for oral cancers in populations of European descent. In this study, we performed the first GWAS of oral cancer in Taiwan with 1529 cases and 44,572 controls. We confirmed two previously reported loci on the 6p21.33 (HLA-B) and 6p21.32 (HLA-DQ gene cluster) loci, highlighting the importance of the human leukocyte antigen and, hence, the immunologic mechanisms in oral carcinogenesis. The TERT-CLMPT1L locus on 5p15.33, the 4q23 ADH1B locus, and the LAMC3 locus on 9q34.12 were also consistent in the Taiwanese. We found two new independent loci on 6p21.32, rs401775 in SKIV2L gene and rs9267798 in TNXB gene. We also found two suggestive novel Taiwanese-specific loci near the TPRS1 gene on 8q23.3 and in the TMED3 gene on 15q25.1. This study identified both common and unique oral cancer susceptibility loci in the Taiwanese as compared to populations of European descent and shed significant light on the etiology of oral cancer in Taiwan.

Cancer carrier screening in the general population using whole-genome sequencing.

BACKGROUND: Cancer is a major cause of death, and its early identification and intervention have potential for clinical actionability and benefits for human health. The studies using whole-genome sequencing (WGS) and large samples analysis of cancer-related genes have been rarely done. METHODS: We performed WGS to explore germline mutations in coding and non-coding areas of cancer-related genes and non-coding driver genes and regulatory areas. Structural variants (SVs) was also analyzed. We used several tools and a subgrouping method to analyze the variants in 1491 healthy participants. Moreover, 275 cancer-related genes sequencing was carried out in 125 cancer patients. RESULTS: The incidence of familial cancer in the Taiwanese general population is 8.79% (131/1491). Cancer carrier rate of cancer-related genes is about 7.04% (105/1491) for pathogenic/likely pathogenic variants (P/LP) on ClinVar database only, and 28.24% (421/1491) for P/LP and loss of function variants. The carrier frequencies of cancer-related genes P/LP on ClinVar database were as follows: 8.40% (11/131), 7.11% (28/394), and 6.83% (66/966) in FC, 1MC, and nMC, respectively. The SVs and non-coding driver gene variants are uncommon. There are 1.54% (23/1491) of actionable cancer genes in American College of Medical Genetics and Genomics (ACMG), and the germline mutation rate of 275 cancer-related genes is 7.2% (9/125) in cancer patients including 4.0% (5/125) of actionable cancer genes in ACMG. After analyzing the frequencies of P/LP variants on GJB2 and SLC25A13 genes, we suggest that these two genes may not be cancer-related genes and need be re-evaluated. CONCLUSIONS: WGS analysis can completely detect germline mutations in cancer carriers. This study use subgrouping approach for samples provides a strategy to study whether a gene or variant is a cancer-related gene or variant in the future studies.

Genetic contributions to female gout and hyperuricaemia using genome-wide association study and polygenic risk score analyses.

OBJECTIVES: To identify genetic variants and polygenic risk score (PRS) relating to female gout and asymptomatic hyperuricaemia (AH) in a genome-wide association study (GWAS). METHODS: Gout, AH and normouricemia controls were included from Taiwan biobank and China Medical University Hospital. All participants were divided into discovery and replication cohorts for GWAS. PRS was estimated according to whether the variant exhibited a protective effect on the phenotypes or not. Each cohort was separated into two groups by the age of 50 years old. RESULTS: A total of 59 472 females were enrolled, and gout and AH occupied 1.60% and 19.59%, respectively. Six variants located in genes SLC2A9, C5orf22, CNTNAP2 and GLRX5 were significantly predictors of female gout in those aged ≥50. For those aged <50 years old, only the variant rs147750368 (SPANXN1) on chromosome X was found. Most variants located in genes SLC2A9, ZNF518B, PKD2 and ABCG2 were found to be significantly related to AH in both age groups. The PRS could explain ∼0.59% to 0.89% of variance of gout in variants with protective effects, which showed 6.2 times of mean PRS in the risk variants, but only 1.2 times in the AH phenotype. Moreover, the PRS also revealed a dose-response trend between AH rates and quartile scores. CONCLUSION: The variants in gene SLC2A9 are the major genetic factors for females associated with gout in those aged ≥50. PRS can provide a more robust prediction of the gout/AH under a homogeneous selection of variants that show effects on the traits.

Polygenic risk score trend and new variants on chromosome 1 are associated with male gout in genome-wide association study.

BACKGROUND: Gout is a highly hereditary disease, but not all those carrying well-known risk variants have developing gout attack even in hyperuricemia status. We performed a genome-wide association study (GWAS) and polygenic risk score (PRS) analysis to illustrate the new genetic architectures of gout and asymptomatic hyperuricemia (AH). METHODS: GWAS was performed to identify variants associated with gout/AH compared with normouricemia. The participants were males, enrolled from the Taiwan Biobank and China Medical University, and divided into discovery (n=39,594) and replication (n=891) cohorts for GWAS. For PRS analysis, the discovery cohort was grouped as base (n=21,814) and target (n=17,780) cohorts, and the score was estimated by grouping the polymorphisms into protective or not for the phenotypes in the base cohort. RESULTS: The genes ABCG2 and SLC2A9 were found as the major genetic factors governing gouty and AH, and even in those carrying the rs2231142 (ABCG2) wild-genotype. Surprisingly, variants on chromosome 1, such as rs7546668 (DNAJC16), rs10927807 (AGMAT), rs9286836 (NUDT17), rs4971100 (TRIM46), rs4072037 (MUC1), and rs2974935 (MTX1), showed significant associations with gout in both discovery and replication cohorts (all p-values < 1e-8). Concerning the PRS, the rates of gout and AH increased with increased quartile PRS in those SNPs having risk effects on the phenotypes; on the contrary, gout/AH rates decreased with increased quartile PRS in those protective SNPs. CONCLUSIONS: We found new variants on chromosome 1 significantly relating to gout, and PRS predicts the risk of developing gout/AH more robustly based on the SNPs' effect types on the trait.

Comprehensive characterization of pharmacogenes in a Taiwanese Han population.

Pharmacogenetic (PGx) testing has not been well adopted in current clinical practice. The phenotypic distribution of clinically relevant pharmacogenes remains to be fully characterized in large population cohorts. In addition, no study has explored actionable PGx alleles in the East Asian population at a large scale. This study comprehensively analyzed 14 actionable pharmacogene diplotypes and phenotypes in 172,854 Taiwanese Han individuals by using their genotype data. Furthermore, we analyzed data from electronic medical records to investigate the effect of the actionable phenotypes on the individuals. The PGx phenotype frequencies were comparable between our cohort and the East Asian population. Overall, 99.9% of the individuals harbored at least one actionable PGx phenotype, and 29% of them have been prescribed a drug to which they may exhibit an atypical response. Our findings can facilitate the clinical application of PGx testing and the optimization of treatment and dosage individually.

Exploring RNA modifications, editing, and splicing changes in hyperuricemia and gout.

Hyperuricemia and gout are two of the most common metabolic disorders worldwide; their incidence is increasing with changes in lifestyle, and they are correlated with many diseases, including renal and cardiovascular diseases. The majority of studies on hyperuricemia and gout have focused on the discovery of the associated genes and their functions and on the roles of monocytes and neutrophils in the development of gout. Virtually no studies investigating the epigenomics of gout disease or exploring the clinical significance of such research have been conducted. In this study, we observed that the expression of enzymes involved in RNA modifications or RNA editing was affected in uric acid (UA)- or monosodium urate (MSU)-treated cell lines. RNA alternative splicing and splicing factors were also affected by UA or MSU treatment. We used transcriptome sequencing to analyze genome-wide RNA splicing and RNA editing and found significant changes in RNA splicing and RNA editing in MSU- or UA-treated THP-1 and HEK293 cells. We further found significant changes of RNA modifications, editing, and splicing in patients with gout. The data indicate that RNA modifications, editing, and splicing play roles in gout. The findings of this study may help to understand the mechanism of RNA splicing and modifications in gout, facilitating the development of new diagnostic and therapeutic strategies.